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I . Abstract. Vacuum-energy calculations with ideal reflecting boundaries are plagued 

by boundary divergences, which presumably correspond to real (but finite) physical 
effects occurring near the boundary. Our working hypothesis is that the stress tensor 
for idealized boundary conditions with some finite cutoff should be a reasonable ad 
hoc model for the true situation. The theory will have a sensible renormalized limit 
when the cutoff is taken away; this requires making sense of the Einstein equation 
with a distributional source. Calculations with the standard ultraviolet cutoff reveal an 
inconsistency between energy and pressure similar to the one that arises in noncovariant 
regularizations of cosmological vacuum energy. The problem disappears, however, 
if the cutoff is a spatial point separation in a "neutral" direction parallel to the 
boundary. Here we demonstrate these claims in detail, first for a single flat reflecting 
wall intersected by a test boundary, then more rigorously for a region of finite cross 
section surrounded by four reflecting walls. We also show how the moment-expansion 
theorem can be applied to the distributional limits of the source and the solution of the 
Einstein equation, resulting in a mathematically consistent differential equation where 
cutoff-dependent coefficients have been identified as renormalizations of properties of 
the boundary. A number of issues surrounding the interpretation of these results are 
aired. 
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1. Energy density and pressure near a plane boundary 

Introductory remarks. The effects of reflecting boundaries in quantum field theory 
continue to be of intense interest, from the mathematical to the experimental [TJ |2j [3j H] . 

From a mathematical point of view there are two kinds of effects. First, there is the 
classic Casimir energy, which is distributed throughout space, associated semiclassically 
with periodic classical orbits, and finite after subtraction of the zero-point energy of 
infinite space. (Here we consider only space-time that is locally flat, and without any 
external potentials that would create additional distributed divergences.) Second, there 
is energy that is concentrated near the boundaries, associated with short paths, closed 
but not periodic, that bounce off the boundaries, and usually divergent in this sense: 
For a perfectly reflecting boundary (such as a perfect conductor in electromagnetism or 
a Dirichlet or Neumann boundary for a scalar field), the vacuum expectation value of 
the energy density has a nonintegrable behavior, 

+ 3 + a) 

s 4 s A 

where s is the distance from the boundary. (The leading exponent is d + 1 when the 
spatial dimension is d.) 

In the original Casimir scenario of the electromagnetic field between parallel flat 
conducting plates [5j [6], energy of the second type is absent. In the famous Boyer 
calculation [7] for the electromagnetic field inside and near a thin spherical shell of 
conductor, the total energy again comes out finite, though only because of a sequence of 
special cancellations [8] (between electric and magnetic modes for c\ , between interior 
and exterior for c<i and C4 , and for a variety of complicated geometrical reasons [HI [10] for 
C3). In general the divergent boundary contributions are present and are playing a major 
role in the ongoing overhaul of the theoretical ideas forced by the huge improvements 
in experimental results since 1997. 

"Analytic" regularization methods (dimensional regularization and zeta functions) 
usually yield finite values for total energy, and ultraviolet-cutoff calculations can reach 
the same values by discarding the divergent leading terms in an expansion in powers of 
the cutoff parameter; but much current opinion tends to limit the scientific validity of 
such methods to calculations of the forces between disjoint, rigid bodies. In studying 
deformable bodies or gravitational effects, the actual energy density near the boundary 
must be taken seriously, together with the energy of the boundary material itself, 
from which it cannot be cleanly separated. It is agreed that divergences result from 
overidealization of the physics of the boundary. For example, a real electromagnetic 
conductor is not perfectly conducting at arbitrarily high frequencies, and at truly 
small scales its atomic structure becomes relevant. For a realistic boundary the effects 
represented crudely by inverse powers in ([1]) are finite, but perhaps large. 

In our present program we attempt to salvage some of the mathematical advantages 
of the old theory by hypothesizing that the vacuum stress tensor, T^ v {x\ of a quantum 
field interacting with an idealized boundary with an exponential ultraviolet or similar 
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cutoff kept finite is a reasonable ad hoc model for a realistic physical situation. We are 
particularly interested in the interpretation of the boundary divergences in the context 
of general relativity, an issue raised in the seminal paper of Deutsch and Candelas [8]. 
Since the actual gravitational effects of vacuum energy in the laboratory are surely very 
small, the problem is primarly one of establishing consistency of the properly interpreted 
theory, so this approach should be sufficient — certainly much more satisfactory than 
simply discarding divergent quantities that are integrals of perfectly finite and plausible 
local densities. The expectation is that the theory will have a sensible renormalized limit 
when the cutoff is taken away; this requires making sense of the Einstein equation with a 
distributional source, interpreting the cutoff as a "regularization" of distributions in the 
mathematical sense [11]. (See also the discussions of mass renormalization in [T2| IT3"].) 

This program was initiated in [14]. Detailed calculations of the stress tensor in a 
model with two spatial dimensions were presented in [15] . where complete consistency 
was claimed between force calculations based on the pressure and those based on 
differentiation of total energy with respect to a geometrical parameter. However, close 
inspection later revealed a sign discrepancy for one pair of terms, those related to 
the leading term in the cutoff boundary energy near a boundary; in dimension 3 this 
problem persists, compounded by a discrepant factor of 2 [16J. There is no evidence 
that this anomaly affects the finite terms constituting the standard Casimir force. 
Nevertheless, since both methods of computing forces are widely used and assumed 
to be equivalent [17] , understanding the phenomenon is imperative. The calculations in 
[15] and [IE] used a exponential ultraviolet cutoff, roughly equivalent to "point-splitting" 
regularization in the time direction. It is reasonable to think that the trouble is related to 
the resulting violation of relativistic invariance. Therefore, a subprogram was launched 
to consider point-splitting in other directions [18]. For rectangular configurations it has 
been found that a cutoff in a "neutral" direction (neither the time, nor the direction 
whose pressure is being sought) yields consistent and physically plausible results. The 
purpose of the present article is to present detailed results along this line, together with 
recapitulation and updating of the material previously published only in conference 
proceedings [HI [16j [18] . 

In this section we review the basic theory of the vacuum stress tensor of a scalar 
field and its coupling to the gravitational field in linear approximation, and we study 
the energy-pressure relations and distributional limits associated to a plane Dirichlet 
boundary. The physical and philosophical issues raised are discussed in section El The 
final section takes a more complete look at the energy-pressure relations in a rectangular 
box (more precisely, a waveguide) when the roles of all four boundaries are carefully 
taken into account. 

Basic equations: quantum field theory. Generalized to curved space-time, the action, 
Lagrangian function, and stress tensor of a scalar field in a cavity Q are 

S = I L ^gd d+1 x, L = — hg^d^dA + m\ T^ = ^-^-. (2) 
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We adopt the metric convention in which g 00 < 0- The parameter £ labels different 
gravitational couplings. In the flat-space limit the field equation and (classical) total 
energy are independent of £, but the stress tensors are different. More specifically, that 
field equation is 



dt 2 



(3) 



accompanied by boundary conditions on Q, which we shall take to be Dirichlet here 
(0 = 0). The simplest formulas for the stress tensor are obtained when £ = |, which 
corresponds to replacing the term (V0) 2 in the energy density for minimal coupling 
(£ = 0) by — 0V 2 0, to which it is related by integration by parts: 

I r/<9</>\ 2 



Too Q) 



Tjj (4) 2 



dt) 



0V 2 <? 



dx 2 . 



for j ^ 0. 



For other values of £ one has 

T^(£) = T^(i)+AT M 

AT 00 = -2 (£-I)V-(0V0). 



- fiV 1 



AT, 



,1.1 



Kg) -E(£ 



dxi 



d 2 



dx 2 



(4) 
(5) 

(6) 
(7) 

(8) 



Henceforth we consider £ = | unless otherwise indicated. The off-diagonal components 
of the tensor play no role in this paper. 

The vacuum expectation values of components of the stress tensor can be expressed 
in terms of the Green function (called cylinder kernel or Poisson kernel) 



00 1 

T(t,r,r') = -£-i 



(9) 



n=l 



where n and u n are the eigenfunctions and eigenfrequencies of the cavity. The formulas 



arc 



P = (Too) 



Pj = ( T jj) 



1 8 2 T 
2~cW 



d 2 T d 2 T 
+ 



dxj 2 



dx'* 



8 2 T 
dxj dx'j 



(10) 



(11) 



where it is understood, formally, that r' is set equal to r and t to at the end. Henceforth 
the angular brackets will usually be omitted from expectation values. 
For more details of the foregoing formalism see |15j . 

In fl9|)- f|TT]) t is an ultraviolet cutoff parameter, not the physical time. However, 
it can be thought of as arising by a Wick rotation from the difference of two physical 

— i(x° — x'°). In the standard ultraviolet-cutoff approach, t 



time coordinates: t 
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is kept different from until the latest possible moment, so that (ITUj) -( TTTj) define a 
nonsingular function even when r' = r. On the other hand, if one thinks of x° as 
just another coordinate in parallel to the components (xj) of r, then it is natural to 
consider implementing a cutoff by separating the primed and unprimed points in some 
other direction in space-time. At the spectral level this maneuver does not lead to a 
classically convergent summation like (Q, and if the separation is in a spatial direction, 
this technical problem persists even if the spatial coordinate in question is subjected 
to a Wick-like analytic continuation. The series will converge in a distributional sense, 
however. Moreover, when a closed-form expression for T is available, as it is in all cases 
studied in this paper, the prescription can be applied to it without encountering any 
divergence. 

In its original application to a theory in curved space-time without boundaries [19] , 
the point separation could be in an arbitrary space-time direction, parametrized by a 
tangent vector tu^. In the particular case of flat space-time the resulting formula of 
Christensen is 

^ = ^f^-4^Y (12) 



2vr 2 t 4 UpU P 

The usual interpretation is that the direction-dependent term averages to 0, while the 
scalar cutoff parameter t is arbitrary; thus the stress tensor of flat space is 

TP = A<T, (13) 

where A is a cosmological constant to be determined by observation. (For an up-to-date 
discussion of this situation see [20] .) Our present concern, however, is with a static 
boundary, and it seems that in this context the only sensible splitting directions are 
pure time (real or imaginary) and all the spatial directions parallel to the boundary 
surface U 

A reflecting wall intersected by a test wall. Consider a plane Dirichlet boundary at 
z — 0. We write rj_ = (x,y) = (xi,x 2 ) for the coordinates parallel to the wall. The 
cylinder kernel for empty Minkowski space is 

T ° = - A T^rr — 7v2 ■ ( 14 ) 

2n 2 t 2 + |r — r'| 2 

The kernel for the half-space z > with the reflecting boundary is instantly obtained 
by the method of images, by adding to T the reflected term 

T = 2^2 t 2 + |r ± - r'J 2 + (z + z') 2 ' (15) 

For determining the physically relevant stress tensor one needs only the term (Tl5|) . 
because the contribution of T is just ( 1T3"|) . ubiquitous and already considered. 



| To foreshadow: From (fT6|) -(|T8 |) it can be seen that if the separation mixes space and real time, 
the denominators will vanish somewhere in the physical region, which defeats the purpose of a cutoff. 
Splitting involving z — z' is ambiguous at best and problematical to define when the point is closer to 
the boundary than the cutoff distance. 
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Figure 1. The reflecting plate is in the x-y plane. The test wall is momentarily in 
the z-y plane. The force on the test wall from the vacuum energy to its left is to be 
studied. See section [2] for discussion of the ignoring of the right side of the test wall, 
and see section [3] for discussion of the neglect of the finite size and other walls of the 
cavity on the left. 



Therefore, we now reserve the notation T for the reflection term alone (sometimes called 
the "renormalized" cylinder kernel). 

We can now calculate the energy density and pressure from ( fTOl) and (TTTj) . After 
the differentiation, without loss of generality we set r' ± = (because of translation 
invariance) and z' = z (because we will not consider point-splitting perpendicular to 
the boundary). Then t and the components of r_i_ are still available as cutoff parameters. 
We define 

M = t 2 + x 2 + y 2 + Az 2 . (16) 
Then the results are 

2ix 2 p = M~ 3 [-3t 2 + x 2 + y 2 + Az\ (17) 

2n 2 Pl = M- 3 [-t 2 + 3x 2 - y 2 - 4z 2 }, 2vr 2 p 2 = M~ 3 [-t 2 - x 2 + 3y 2 - 4z 2 }, (18) 

and P3 = 0. Qualitatively, one can say that in the vacuum expectation there is a layer 
of energy against the wall, where M is smallest. The vanishing of p% is understandable, 
because a rigid perpendicular displacement of the wall does not change the total energy 
of the configuration. 

To understand the nonvanishing pressures parallel to the wall, imagine another 
planar boundary — a "test wall" — at x = and find the force on it (from the left side 
only — see figured]). The volume of space occupied by boundary energy increases with 
x, so the total energy on the left also changes linearly when the test wall moves. 

In accordance with the principle of energy balance (also called principle of virtual 
work) one expects 

dU f°° 

Pi dz = — — — — — I pdz = —E. (19) 
o ux J Q 

Here a trivial transverse dimension is suppressed in the notation: F is a force per unit 

length, E is an energy per unit area, and U is the total energy per unit length, equal 

to E times the length of the cavity (which may be arbitrarily large). If all the cutoffs 

are removed, one has p = (32it 2 z 4 )~ 1 = —p 1 , so energy balance is formally satisfied 

point-by-point in z, but of course the integrals are divergent. 
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Consider now the traditional ultraviolet cutoff, t ^ but r_|_ = 0. Evaluating the 
integrals in f|T9|) with the integrands from (fTTl) and ([TBI , one finds that 

a clear discrepancy with f|T9|) in both sign and magnitude. (For details see the 
Appendix.) The E defined in this way is the same as obtained from the second (surface) 
term in the small-t expansion of 

e = - li r ' r) d " x = \ % L ° ne ~ tujn - (2i) 

But we shall argue that this E is wrong, whereas the F and (1191) are correct (within the 
framework of the model). 

Next look at point-splitting perpendicular to the test wall (x ^ 0, t = = y). It is 
easy to see that the calculation is identical to the previous one except that (t, p) change 
places with (x, —pi). Thus the result is 

F = +2E. (22) 

Finally, consider point-splitting in the "neutral" direction, y ^ 0, t = = x. One 

gets 

F = —E, (23) 

as should happen, according to (Tl9l) . In fact, the balance relation holds pointwise for 
the integrands; their values are 

9 = 2n 2 (y 2 + Az 2 ) 2 ' Pl = ~ 2<K 2 (y 2 + 4z 2 ) 2 ' ^ 

The quantities in (j23|) and fpMl) agree with (|20|) for the pressure and with (!22l for the 
energy. We therefore propose to adopt them as "correct". Note that for the Dirichlet 
wall, this E is positive, unlike the traditional boundary energy from f l2T]) . 

A similar calculation for the "conformal correction" terms (J7j)-(JED is entirely 
uneventful. They do not exhibit the anomaly: Api = —Ap always. Furthermore, 
as long as the cutoff is finite, their integrals over all z are anyway, as befits their 
origin in an integration by parts in a theory where the field vanishes on the boundary. 
(That this consistency relation fails without the cutoff was pointed out by Ford and 
Svaiter [21].) 

It may appear strange that the pressure and the energy have opposite sign. 
This behavior is standard, however, for an energy density that remains constant as 
a geometrical parameter is varied, so that the total energy is proportional to the volume 
occupied. It holds, for instance, for the surface tension between two fluid media [22j 
pp. 84-87], and also for cosmological "dark energy" (cf. (ITBl). 
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Basic equations: general relativity. In the notation of [23], the linearized Einstein 
equation is 

□V = , (25) 

where 

hp, = V/ - — (Tr h)77^ , h^ = - rj^ . (26) 

We apply (|25|) to our model of a quantized scalar field with a single plane boundary, the 
boundary condition being one of the point-splitting cutoffs discussed above, approaching 
the Dirichlet condition as a limit. Gravitational effects in the lab, although formally 
infinite in the theory without cutoff, are presumably actually tiny. Therefore, the 
linearized Einstein equation (with a flat background) should be physically adequate. 
In a slightly better approximation one might use a curved background to represent the 
effect of the mass of the boundary. 

We assume a static solution. Then, tentatively adopting the usual ultraviolet cutoff, 
we have the basic equation of [Tj 



8 4r 2 _ 3f 2 

- V 2 /,oo = 16^ = - {t2 + Az2f 9{z) + (■ ■ .)6{-z). (27) 

Because we consider only the region z > and the equation is linear, the term describing 
matter sources at negative z decouples and can be ignored. However, it is essential to 
include the factor d{z) to study the singular limiting behavior of /ioo and p at z = 0. 
For an infinite plane wall, h 00 should depend only on z, so V 2 h 00 = h,Q (z). 
The solution of (1271) that vanishes for z < is 

4z _i/2z\ 1 1- 



h 00 (z) - 9{Z) 



, , (28) 

For a complete solution one should add a similar term with support in z < representing 
the effect of the matter in that region, and also a possible global solution (affine linear 
in z) of the homogeneous equation. (If the coefficient of z in this homogeneous solution 
is made different on the two sides, a delta function representing the gravitational field 
of the wall itself will be added to p.) 

Taking the limit t 4- in (1271) will yield a differential equation with a distribution as 
its source. Taking that limit in the solution will yield a singular distribution. Both limits 
involve somewhat arbitrary regularizations (Hadamard finite parts). Before presenting 
the details, we review some needed distribution theory in the next subsection. 



Distributions and the moment expansion theorem. Our analysis requires the study of 
the behavior of several quantities of physical interest as a parameter, t, approaches 0. 
The singularity of the quantity suggests that instead of an ordinary limit one needs 
to consider a distributional limit. Fortunately, a distributional theory of asymptotic 
expansions is available |24|. 

The key result in the distributional theory of asymptotic expansions is the moment 
asymptotic expansion, that says that if / is a distribution of one variable, / e V (R) , 
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that decays rapidly at ±00, then the asymptotic behavior of / (Ax) as the parameter A 
becomes large is 

n=0 

Here the {/i„} are the moments of the distribution /, namely, 

/i n = (f(x),x n ) , (30) 

which become fi n = f (x) x n dx if / is a locally integrable function or become 
= J2T=-oo a k(bk) n if / (x) = J2T=-oo a k5(x - b k ) is a train of delta functions. 

Observe that the behavior of / at ±00 is given by considering its parametric 
behavior; this reflects the fact that distributions do not have point values, in general, 
as ordinary functions do. An alternative approach that uses the ideas of Cesaro 
summability is possible [25], but in the present article it is exactly the parametric 
behavior that we need to consider. 

We used the term "of rapid decay at ±00" to describe the distributions for which 
the moment asymptotic expansion holds. In technical terms, / is of rapid decay at ±00 
if it belongs to the distribution space JC' (R) , dual of JC (R) . The space K is formed by 
the so-called GLS symbols [26J; a smooth function <p belongs to JC if there is a constant 
7 such that (f)( k \x) = 0(|x| 7_fc ) as |x| — > 00 for k = 0, 1, 2, . . . , that is, if 0(x) = 0(|x| 7 ) 
strongly; the topology of JC is given by the canonical seminorms. 

Notice that if £' (R) is the space of distributions with compact support, then 
£' (R) C JC' (R) , so that the moment asymptotic expansion holds for distributions of 
compact support. 

If / is a distribution of one variable that does not belong to JC' (R) then it does 
not satisfy the moment asymptotic expansion, and, in general, the behavior of / (Ax) as 
A — > 00 is very complicated. There is one situation, however, when / (Ax) has a simple 
expansion [27] . Indeed, suppose that / has support bounded on the left, supp / C [a, 00) 
for some a, and suppose that / is an ordinary locally integrable function for x large and 
that the ordinary expansion 

/ (x) = b x x^ + • • ■ + b n x Pn + O (x p ) , asx-^oo, (31) 

holds, where f3\ > (3 2 > ■ ■ ■ > f3 n > (3 and where — (k + 1) > j3 > — (k + 2) for some 
integer k. Then / (Ax) has the distributional development 

/ (Ax) = ± (Ax) + £ ■, ffif {X) + O (AO , (32) 
j=l j=o J ' 

as A — > 00. Here the distributions gj (x) are suitable regularizations of the functions 
9 (x) x^ j , namely Vf (0 (x) x^ ) , as we explain below, while the /ij are still the moments, 
but understood in a generalized sense, since the integrals J a °° / (x) x^dx will be divergent 
at infinity, in general, and thus we need to consider their Hadamard finite part, 
F.P. f a °° / (x) xHx. 
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Interestingly, the asymptotic expansion (1321) considers the behavior of / at infinity, 
but it is writen in terms of distributions that could be singular at the origin. The process 
of associating a distribution to a singular, not locally integrable, function is considered 
in [11] ; in the mathematical literature one says that the distribution is a regularization 
of the singular function. There are many regularization procedures, and, in general, 
regularizations are not unique |28j; for our purposes we need to consider the Hadamard 
finite part of integrals and the corresponding regularization method. 

Let us first recall that the finite part of the limit of F (e) as e — > + exists and 
equals A if we can write F (e) = F fin (e) + F infin (e) , where the infinite part -F infin (e) is 
a linear combination of functions of the type e~ p In 9 e, where p > or p = and q > 0, 
and where the finite part Fg n (e) is a function whose limit as e — > + is A [23]. We then 
write 

F.P. lim F(e) = A. (33) 

£^0+ 

A similar idea can be applied to limits at infinity, which turns out to be equivalent to a 
change of variables, 

F.P. lim G (A) = F.P. lim G (1/e) . (34) 

Suppose now that we have an integral J f (x) dx such that J f (x) dx is well defined 
for c > a, but for which f^f(x)dx is divergent. Then one considers the function 
F (e) = J^ +£ f (x) dx for e > 0, and if the finite part of the limit exists as e — > + , then 
we call this value the Hadamard finite part of the integral, and write 

rb rb 

F.P. / / (x) dx = F.P. lim / / (x) dx . (35) 

J a e ^ 0+ Ja+e 

Let / now be a function of one variable that vanishes for x < a, that is locally 
integrable in the open interval (a, oo) , but for which the integral f (x) dx is divergent 
for b > a. Then, in general, one cannot associate a distribution in T>' (R) to /. 
Nevertheless, if the Hadamard finite part of the integral, F.P. J f (x) dx, exists for 
b > a, then one can define a distribution Vf (/) , its "partie finie", by putting 

(Vf(f(x)),<P{x)) = F.P. I f(x)<P(x)dx, 0GX?(R). (36) 

J a 

In the case when f (x) =6 (x) x 13 for some /3 G C we obtain the distribution 
Vf(9(x)xP). When 3?e (3 > —1, however, the function 9 (x) x 13 is locally integrable 
even at x = and thus Vf {9 (x) x 13 ) reduces to the standard distribution associated 
to 9 (x) x 13 ', which is usually denoted by x + \ one can perform analytic continuation of 
x+, and obtain a distribution for (3^—1, —2, —3, . . . , and it can be shown [24] that 
Vf {9 (x) x 13 ) = for such values of (3. On the other hand, the symbol x + does not make 
sense for (3 = —1, —2, —3, . . . , while Vf (9 (x) x~ k ) is well defined for k = 1, 2, 3, 

The formulas 

x()' = f3xl~\ (37) 
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(Ax)J = AVJ , (38) 

hold for p ^ — 1, —2, —3, . . . since they hold for $fte/3 large, by the principle of analytic 
continuation. However, we have the modified formulas 

'O(x) 



dx 


6(x] 


d 

dx 


( vf 



X 



(39) 



6(x)\\ ,„„fe(x)\ (-l) k S^ 



kVf[^) + " ~' kX) , (40) 



and 



x k ) J V x k+1 J kl 



«|W -I w ffl] + H^W (41) 



Distributional limits of energy and pressure. In [M] the foregoing mathematical theory 
was applied to fl27D and Let X — l/t. The limit of (J27J) is 



dx 2 2vr 



Z^ 71 



+ sl r(2) -I^ ln(2A)r(2) + (l)- (42) 



and the limit of (128]) is 



/loo(^) = 2A 3 fl(^ - -A 2 0(z) - — F.P. 



7T 127T V, Z 2 

(The apparent dimensional incoherence of the 6' ln(2A) term is attributable to the 
violation of scale invariance by the finite-part operation.) 

One quickly verifies that fj4*3|) is a solution of ( 142]) . The limit of the solution is 
the solution of the limit. This also follows from the definitions of the distributional 
operations. 

A consequence of the moment expansion theorem is that, when one uses a 
singular function to define a distribution, divergent leading powers can be replaced 
by derivatives of 5 with arbitrary finite coefficients [IT]. These terms can be interpreted 
as renormalizations of properties of the boundary. The previous observation thus shows 
that our toy Einstein equation survives the renormalization process as a mathematically 
consistent differential equation. 

We shall now go beyond [H] in two ways. First, keeping t as the cutoff parameter, 
we examine the differential equation for the pressure component, using fflBl with 
x = = y: 

d 2 h u (z) ^ 8 1 „. . , AA . 

-^ i = 16 ^ = -?(FTW eW ' (44) 

The appropriate solution analogous to ff28|) is 



hn(z) -- 



IT It 3 



2z . _ x f2z 



tan 



(45) 
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For the evaluation of the asymptotic behavior of (1441) . the relevant distribution to 
use in the moment expansion theorem is 

^(irW^ii^ ^)^ 00 ' (46) 

The moment expansion theorem states, up to the relevant order, that the asymptotic 
expansion of fi(\z) is 

(47) 



16 1 A 4 ' V z" I 3!A 4 

3=0 J ' 

where the moments ^j(fi) of the function f\ are 

Mfl) = I M I \ dz = I > ( 4J 



1 


(1- 


f 4z 2)2 




1 


(1- 


f 4z 2 ) 2 




1 



M/i)= / n I 7I^ zdz = -, (49) 

^ (/i) = / (irir" 2 ^ = ^' (50) 
f°° i ii 

"^)=F. P .y o irTW ^ = __ + _ ln2 . (si) 

Therefore, the distributional limit of the differential equation (|44p is 

rf 2 /iii(2;) 



1 5"'(z) + — \n(2\)8"'(z) + o(\). (52) 



dz 2 



24vr v ' 12tx y ' y ' \\ 
Similarly, the relevant function for the analysis of ( 143]) is 

7T 1 11 „ / 1 



/ 2 (s) = ztan- 1 (2z)^(x) = |* -- + _- + C>(-). (53) 
The moment expansion theorem says that 

/2 (A,)~^(A,)(Az)-i^) + ip/(H 



7! 



^^-^w + ^(^)-^ kW 'w 



i=o 
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Finally, the relevant moment expansion coefficients this time are 

f°° 7T 

fMj(h) = / ^tan- 1 (2z) dz = — , (55) 



o 



16 



poo 11 

/ z 2 ttm- l {2z)dz = — + — In 2. (56) 
J 1 ; 72 24 v ; 



Forming the correct linear combination of these terms, we have 

h n (z) = \ 3 6(z)z - -X 2 6(z) + + ^(*) 

7T 1Z7T V Z* / 8 



1 1 

- + — In 2 
L72 24 



*' W -J-.„Ay w+0 (i 



71 YllX \ Z Z / 8 367T 

-±-H2X)5'(z) + oQ. (57) 

By taking the second derivative of the solution (1ST)) , according to the Hadamard 
T 3 / formulas, we verify that the differential equation ( 152")) is satisfied. This completes 
the counterpart for the pressure of the energy calculation in [T4] . 

Next we want to see what happens when we use y as the cutoff parameter. It is 
neutral with respect to both the time-energy and the x-pi dimensions. The pertinent 
formulas are (I24p . One immediately sees that the pi calculation is identical to what we 
just did, with y in the role of t, and also that p in this case is just the negative of pi . 
So we do not need to do any more calculating, and we have three new, agreeing results 
(two for pi and one for p) that outvote the old formula for p in [H] and l )2TI) . Similarly, 
we have formulas for Tin and hoo that now seem more trustworthy than the old formula 
for h 00 in HM and (1281). 



Comparison of old and new energy formulas. The old formula (I43p for the distributional 
limit of hoo, based on the t cutoff, was 



h° d (z) = 2\*e(z)z - -x'e(z) - ^-vf(- 

n YLix \ 



z 2 



' 5'{z) + —\tl{2X)5'{z). (58) 



187T ' ' 127T 

The new one, based on the y cutoff, is the negative of ((5 



/^) = -\*0(z)z+ h 2 9(z) - ±-Vf{^) - \\5{z) 

1 rJ, x 1 



367T 127T 

The structure is exactly the same (except for the accidental absence of a 5(z) term 
from (158")) ). Only the numerical coefficients are different; their ratios are simple integers, 
shown in table HJ The coefficients are the same for the regularized 8(z)/z 2 term, which 
is the only term describing a spatially extended vacuum effect of the quantum field 
theory. All the other terms (apart from the bare S'(z) terms, which could in each case 
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Table 1. Ratio h\ 



old /L: 
00 / I 



new 



00 



Term 



Ratio of coefficients 



A 3 z(9(z) 
A 2 #(z) 
A(5(z) 
ln(2A)<5'(z) 
8>{z) 



-2 







f 



2 



1 




f 



be absorbed by rescaling the argument of the logarithm) involve powers or logarithm 
of the large cutoff parameter A. The spirit of renormalization is to regard these terms 
as part of the wall, not the field, and to replace the A-dependent factors by arbitrary 
finite constants. Note that the two terms involving 6(z) can be regarded as part of the 
solution of the homogeneous Einstein equation, and the other two terms are localized 
on the wall. A full discussion of the homogeneous solution is beyond the scope of this 
paper. 

2. Caveats and implications: A dialogue 

Salviati: Friends, I am perplexed by the paradoxical behavior of the energy and pressure 
in my model. The change in energy as the side of a box moves does not match the force 
on the side. Even the sign is wrong! 

Simplicio: I don't understand why you are so upset. We know that the divergent 
terms in vacuum energy are unphysical and must be discarded. The divergences you're 
looking at are not logarithmic, so in an analytic regularization scheme (dimensional 
regularization or zeta functions) they do not arise in the first place. Also, those methods 
are Lorentz-covariant, whereas point-splitting in any particular direction is not; it's not 
surprising that the latter gives noncovariant results for the cutoff-dependent terms. 

Sagredo: I disagree with Simplicio about the correct way to think about divergences 
in Casimir calculations. (With a realistic boundary the "divergences" will be finite but 
nonzero, and they have physical content.) But I, too, think that your "paradox" is 
a nonproblem, for a whole slew of reasons. First, in the calculation symbolized by 
figure [T] you have completely ignored the space to the right of the test wall. Surely it is 
obvious that the force F is cancelled by an equal and opposite force from the other side. 
Second, your scalar field is a toy that doesn't exist in the real world. Nature's massless 
field is electromagnetism, and one of the most famous properties of electromagnetic 
Casimir energy is that the surface divergences you're studying vanish there, because the 
electric and magnetic contributions cancel. Also, there is no surface divergence for a 
flat boundary and a scalar field if £ = |. Third, your cutoff is not a physical model of a 
realistic boundary; it is just an ad hoc procedure. If you use a noncovariant or otherwise 
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unphysical cutoff, you must be prepared to insert unphysical counterterms; they are not 
a physical pathology but an artifact of a suboptimal formulation of the problem. 

Salviati: Give me a few minutes to respond to this barrage! First, this anomaly 
was first noticed in calculations for a thin reflecting spherical shell [S. Fulling and M. 
Schaden, unpublished]. In that case the force in question acts on the sphere itself, not 
a test surface perpendicular to it. The leading term in the surface energy per unit area 
is the same for convex and concave surfaces, and the area (hence the volume occupied 
by the surface energy) depends on radius in the same way. So from the energy point 
of view, the inside and outside forces should be in the same direction! And indeed, the 
direct calculation of force from pressure gave that result, but with that same discrepant 
factor, —\ [16J. 

Simplicio: Why not just show us the sphere calculations? 

Salviati: They are rather complicated (and computer-assisted) and not yet available 
in publishable form. Instead, I have shown you that the same phenomenon occurs in the 
much simpler rectangular configuration. My belief is that energy-pressure balance ought 
to hold for each side of the test wall separately, since they are effectively independent 
systems (especially if the test wall is itself another perfectly reflecting wall). I don't 
know or care what exists on the other side of the wall. But if you insist that I analyze 
a thin partition, taking both sides into account, I can imagine scenarios where the two 
sides don't balance. Recall that the surface energy for a Neumann plate has the opposite 
sign from a Dirichlet plate. So if the plate in figure [1] has the Neumann property for 
x > 0, the resulting force on the test wall will be in the same direction as that from the 
Dirichlet plate on the other side. 

Sagredo: Whoa! Does your plate magically change from one material to another as 
the test wall moves? 

Salviati: At the instant portrayed in figure [Q, the force is as I've described. What 
happens after the test wall has moved a finite distance is not really relevant. But 
anyway, I said I can "imagine" these situations, not create them in a lab. Perhaps the 
reflecting media stretch and shrink as the test wall is moved, but in such a way that 
their reflective properties don't change. I see no reason why such a situation could not 
exist in principle; we're talking here about preserving a very general physical principle. 
Or perhaps the Neumann material is a thin strip that winds around a spool as the test 
wall moves to positive x. 

Sagredo: But the boundary of the Neumann substance still exists even if it is wound 
into a spiral. This reminds me of your somewhat ill-fated pistol design [15]. Leaving 
aside the fact that the pistol sucked in its bullet instead of firing it, you had to admit 
that two conductors in close contact are not the same thing as a solid block of conducting 
material; rather, one has a limiting case of a small gap between conductors, and in the 
gap the Casimir energy is very important. 

Salviati: Yes, but our point was that the idealized theory simply must break down 



Surface Vacuum Energy in Cutoff Models 



16 



when the gap is very small (atomic dimensions), even if we don't know exactly how. 
I admit that the spool scenario has complicating factors. How about this one: There 
is no boundary surface at all along the z axis at positive x. Instead, (to use electrical 
terminology) the right side of the apparatus is filled with a liquid conductor — mercury 
- which is pushed out of the way in the z direction into a large tank as the test 
wall moves. Then there is no boundary vacuum energy in the right half. My point is 
that there is a vast variety of things that could exist on the right of the wall, and no 
way to prove that they all give the same force on the wall. Anyway, this rectangular 
scenario is just a simple model of a more general problem that can't be argued away. 
I have already reported our result about the sphere. My reply to the point about 
electromagnetic surface energy is similar: The leading surface terms cancel there, but 
for curved surfaces there are higher-order terms that don't, and the same issue will 
probably arise. 

Simplicio: You said something that bothered me. You want to make your test 
wall into a perfect reflector in its own right. So doesn't its vacuum energy need to be 
counted? Also, you started by calling your apparatus a "box", implying that it has 
walls on all sides instead of extending to infinity as it does in your figure. But then 
your argument for the correctness of the equation p^ = breaks down. You argued that 
moving the boundary at z = would merely push the vacuum energy around instead 
of changing its total value. That isn't true if the cavity's length is changing. And your 
argument for ( |T9l) explicitly assumed that the cavity's length in the x direction is finite 
and changing! You can't have it both ways. 

Salviati: You are right. For a convincing and complete treatment we should consider 
a bounded system — a rectangular box — and study the energy and pressure associated 
with all those walls. In fact, that has been done (see section [3]). The analogue of the p% 
argument does work out, but contributions from infinitely many image sources must be 
combined to get energy-pressure balance. 

Sagredo: How seriously do you take your model physically? 

Salviati: I don't know. There are two philosophies one could take. First, one might 
hope that the finite-cutoff theory (possibly including some finite, physically motivated 
counterterms) could be a physically plausible model of a real boundary, including its 
gravitational effects, at least at distances not extremely close to the wall. The prototype 
I have in mind is the "hard core" in the Lennard- Jones potential, which nuclear 
physicists regard as adequate when they don't need to probe the chromodynamical 
structure of nucleons at short distances. Surely in that situation a theory that violates 
energy-pressure balance is unacceptable. Neutral point-splitting gives a theory that is 
acceptable in this respect, but I can't regard it as a logically sound, long-term solution, 
since it amounts to changing the rules of the game in each situation until I get a result 
that I like. Its only justification is that, unlike less contrived alternatives, it does not 
immediately produce results that are obviously wrong. In the long run, detailed study 
of the physics of real boundaries will be necessary. 
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The other point of view is that the cutoff is a mathematical means to an end, 
which is a limiting theory with the cutoff removed. At the intermediate stage a 
violation of energy-pressure balance may be tolerated, so long as the final theory 
is physically acceptable. As you both remarked, an unphysical regularization can 
necessitate unphysical counterterms, so that the undesirable terms appear in the final 
equations with coefficient (i.e., not at all). I could live with that, but it is certainly 
highly unaesthetic. 

Simplicio: So, you have left us all with several questions to think about. 

• Are both of those philosophies physically tenable? Or one, or neither? 

• What is the meaning of terms proportional to 5 and 5' in a component of a metric 
tensor? Do we learn anything more than if we had just discarded them? 

• When you put in numbers, is the vacuum stress sufficiently small to justify your 
use of linearized Einstein equations? Must you add counterterms to assure that? 

• In particular, I am concerned by the terms in that are linear in z. What 
is the significance of the linear and constant terms (solutions of the homogeneous 
equation)? To what extent can one exclude them by boundary conditions or remove 
them by gauge (coordinate) transformations? 

• Naively I would expect a term that is independent of A to be "physical" and hence 
well defined. But there is one such term in table [D for which the ratio is not unity. 
Should we be worried about that? I can see that that term could be regarded as 
part of the logarithmic term, although the coefficient ratio for the logarithmic term 
itself is unity. 

• What will you do near the corner of a parallelepiped, where there is no neutral 
direction? 

3. Energy density and pressure inside a rectangular box 

In this section we give a more careful accounting of the total energy associated with 
two perpendicular boundaries as in figure [TJ The test wall itself is made into a perfectly 
reflecting boundary, and to make all integrals proper two other such walls are introduced 
parallel to the original ones. The cavity is thus a finite box as shown in figure [2] More 
precisely, the cavity is an infinite rectangular waveguide, because we have no need to 
introduce boundaries parallel to the x—z plane. Therefore, in the following, "energy" 
should be read as "energy per unit length" (in the y direction), and similarly for "force". 

The wave equation in this space can be solved by a sum over infinitely many images. 
The analysis follows rather closely the two-dimensional theory in [151 E] , but the details 
of the expressions are different. As shown in figure [21 there are four classes of images. 
Points periodically displaced (with periods twice the length, a, and width, b, of the box) 
from the source point at r' = (x', 0, z') contribute to the cylinder kernel T(t, r, r') the 
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Figure 2. Image points of a point x in a rectangle (cf. [HI [T31 [H]), classified as 
periodic displacements (p), horizontal reflections with (possible) vertical periodic drift 
(h), vertical reflections with horizontal periodic drift (v), or corner reflections (c), 



Tpmn 2vr t 2 + y 2 + (2ma + x- x') 2 + (2nb + z- z') 2 ' ^ 



terms 



The term Tp o is the free kernel T and hence should be omitted as it was in (1151) . Points 
reflected horizontally (i.e., reflected through a vertical line in the lattice of figure |2]), 
together possibly with a periodic displacement in the vertical direction, contribute 

T Hmn = + — t2 + y2+ ( 2mG _ - _ x ,y + ^ + z _ z ,y ■ (61) 

Similarly, vertical reflections produce 

— 1 1 

Vmn - + — ~ 2 ~ y2 ~ + x _ ~ _ z _ ■ (62) 

Finally, reflections in both dimensions give rise to terms 

T Cmn = - — t 2 + y2 + ( 2m a -x- x') 2 + (2nb -z- z') 2 ' ^ 

(Such a point is connected to the source point by a straight line that passes through an 
intersection of the horizontal and vertical lines of the lattice. When folded back into the 
box, this line can be characterized as a classical path or optical ray that is reflected from 
a corner of the box. All the other terms have similar classical-path interpretations.) 

On the plates at x = and x = a each term contributes a pressure in the x direction 
according to (11 ip (with j = 1, which will be henceforth understood). One finds that 
(for any x) 

PRmn = = PCmn, (64) 

8fn 2 a 2 1 
PPmn ~ ~ TT 2 (t 2 + y 2 + Am 2 a 2 + An 2 b 2 ) 3 + 2n 2 (t 2 + y 2 + Am 2 a 2 + An 2 b 2 ) 2 ' { } 

Sm 2 a 2 1 
PYmn = ix 2 {t 2 + y 2 + Am 2 a 2 + 4(z - nb) 2 f ~ 2n 2 (t 2 + y 2 + Am 2 a 2 + 4(z - nb) 2 ) 2 ' (66 ^ 

Note that pp mn and p\ mn are independent of x. 
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Similarly, the energy density due to each term can be computed by (imp . For certain 
purposes, however, it is convenient to take only one t derivative before evaluating one 
of the sums or integrals in closed form. In particular, consider (at x' = x, z' = z) 

rp __ dT}imn t (f,7\ 

Hm „- m - ■ -^2^2+^2 + 4 ( x _ ma )2 + 4n 2 6 2)2- 



One has 







n 2 (t 2 + y 2 + Ax 2 + An 2 b 2 ) 

m=—oo " " - ■■ - — ^ ' ' " " 



7r 2 (t 2 + y 2 + Ax 2 + An 2 b 2 ) 2 ' ^ 
which is independent of a. Therefore, the H-energy in the box is unchanged whenever 
one of the vertical walls moves, which is precisely consistent with from (1M1) . Note, 
however, that it was essential to consider all values of m at once; for m = 0, for 
instance, the energy in the box is not exactly constant because a part of the surface 
energy distribution created by one wall is pushed through the opposite wall. 

The V-energy density is easily seen to be independent of x, since, in analogy to (JHj 



Tvmn n 2 {t 2 +y 2 + 4m 2 a 2 + ±{z-nb) 2 ) 2 ' ^ 
Therefore, the total V-energy is equal to that density times a; so its derivative with 
respect to a has two terms, one just equal to the energy density, 

1 9Ty mn /r 7r .\ 

PVmn = ~ 2 ~^f- , (70) 

and the other equal to 

a J^Pvmn • (71) 
One calculates 

2_e 1 

PVmn ~ K 2 (t 2 + y 2 + Am 2 a 2 + A{z-nb) 2 y + 2<K 2 {t 2 + y 2 + Am 2 a 2 + A{z-nb) 2 ) 2 ' ( ' 
If t = (the neutral cutoff), then (1711) equals the negative of the first term of (1661) . and 
(ITOj) is the negative of the second term. Thus the principle of energy balance (fT9l) holds. 
If t ^ and y = (the ultraviolet cutoff), the first term in ( 1721) survives and creates 
the pressure anomaly (specifically, from the terms with m = 0). 
A calculation similar to ( 1681) yields the formula 
t 2 - v 2 

Stt^ + t/ 2 ) 2 1 J 

for the total C-energy. It is independent of a, hence consistent with (I64p . In dimension 
2 the corresponding term vanished entirely. The difference is that in dimension 3 the 
"corner" is really an edge of a parallelepiped, whose true corners would indeed give no 
energy in a cutoff theory. 

The P terms are nondivergent, and of course for p they give the standard (finite 
and constant) bulk Casimir energy density of the waveguide. 
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Appendix A. The pressure anomaly in dimension d 

When the spatial dimension is d (> 2), the free cylinder kernel is, up to a constant 
factor, 

T oc (t 2 + |r-r'| 2 )" ( ^ 1)/2 . (A.l) 
Thus the reflected kernel generalizing (|T5|) is 

T oc (t 2 + \v ± - r'J 2 + (z + z 'f)-^ d - 1)/2 . (A.2) 
Therefore, following (fT0|)-(fTT|) and (jI5j)-(fTg|), one finds 

p oc -M(d)~ (d+3)/2 (~t 2 d + r ± 2 + Az 2 ) , (A.3) 

Pl oc M(d)~^ 2 {t 2 - X 2 d+J2 4 + 4z 2 ) , (A.4) 

with z = Xd , x = X\ , and 

M{d) = t 2 + r ± 2 + Az 2 . (A.5) 

The missing numerical constants in (1A.3I) and (I A. 4 1) are the same, namely [30] 

'd+l' 



l -c{d) = i7r-( d+1 )/ 2 r (i 



(A.6) 



We are to choose one of the coordinates, t or Xj (j ^ d), to serve as the cutoff and 
to set the others to 0. Call that coordinate w. Then we are confronted with integrals 
of the type 



—uj^d _|_ 4_2 2 

Iw= I ^ ( w 2 + 42 2)(d+3)/2 ( A - 7 ) 



and 

w 2 + 4z 2 '' x 



The form (1A.8[) applies if w is "neutral" with respect to the component of energy or 
pressure concerned; (I A. 71) applies if w is not neutral (w — t for energy density, w = x 
for pi , etc.). 
Note that 

I w = I-(d + l)w 2 ^ dz {w2 - 4z2)(d+3)/2 ■ (A.9) 
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Because 

d z (d + l)w 2 d 

dz («j 2 + 4z 2 )(«)/ 2 = (w 2 + Az 2 )( d + 3 )/ 2 ~ (w 2 + Az 2 )^ 1 )/ 2 ' 
the second term in (1A.9|) equals 

00 d 

dZ ( w 2 + 4jZ 2)(d+l)/2 = ~ d ' L 

Thus 

I w = -(d-l)I. (A.10) 

When d = 3, flA~T0l) yields the anomalies in fl20H and fl22D, and (IA~8l) yields the value 
in ( 1201) . When d = 2, the anomaly reduces to a sign change, which was overlooked in 
the discussion near the top of p. 18 of [15J. 
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